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TWO-DIMENSIONAL  NUMESICAL  SIMJIATION  GT  NEGATIVE 
DIFFE31ENTIAL  MOBILITY  SEMICONDUCTING  DEVICES 


T.  M.  McBigh  and  H.  L.  Grubln 

United  Technologies  Research  Center,  East  Hartford,  Connecticut  06IOB 


ABSTRACT 

A nunerlcal  simulation  In  two  spatial  dimensions  of  nega- 
tive differential  nobility  semiconducting  devices  Is  described. 
The  simulation  models  time  dependent  operation  of  two  or  three 
terminal  devices  In  an  external  circuit. 

The  mathenatlcal  model  of  the  semiconductor  consists  of 
a set  of  nonlinear  second-order  partial  differential  equations 
which  are  solved  on  a rectangular  mesh  using  explicit  finite 
difference  techniques.  The  external  circuit  is  modeled  using 
lUB^d  linear  elements  and  the  resulting  equations  are  solved 
using  Implicit  techniques.  The  external  circuit  equations  are 
coupled  to  the  semiconductor  equations  by  the  current  flow  and 
potential  at  the  device  terminals.  The  solution  of  the  result- 
ing time  dependent  boundary  value  problem  is  a feature  of  the 
simulation  which  permits  more  realistic  numerical  analysis  of 
HDM  devices  than  was  previously  possible. 

The  simulation  is  implemented  using  interactive  computer 
graphics.  As  the  solution  evolves,  displays  of  the  potential, 
charge  density  euid  current  flew  in  the  semiconductor  are  avail- 
able  via  a computer  graphic  terminal.  Sample  results  of  simu- 
latiods  of  a negative  differential  mobility  field  effect  traxw 
•later  are  presented. 

DflRODUCTION 

A field  effect  transistor  (abbreviated  FET)  is  a device  in 
which  the  bulk  conductance  of  «ui  extrinsic  semiconductor  is 
controlled  by  an  externally  applied  bias.  .Several  mechanisms 
may  be  used  to  achieve  this  control ; one  of  which  employs  a 
reverse  bias  Schottky  contact  to  control  the  cross-sectional 
area  through  which  current  flows.  A prlncijial  feature  of  fET 
operation  is  that  it  is  a majority  carrier  effect.  The  minor- 
ity carriers  do  not  effectively  participate  in  transistor 
action. 

On^  possible  configuration  of  an  FET  is  shewn  in  Fig.  1. 
The  device  consists  of  an  n-ty7>e  semiconductor  with  low  resis- 
tance contacts  at  either  end,  labeJXed  source  and  drain,  on 
the  top  of  the  semiconductor  is  a Schottky  metaJ  electrode, 
called  the  gate  contact,  which  has  the  effect  of  creating  a 
region  that  is  relatively  free  of  mobile  carriers  (shaded  re- 
gion). The  region  beneath  the  depleted  region  contains  the 
mobile  carriers  and  is  called  the  conducting  channel.  The 
device  generally  sits  atop  of  a semi- insulating  substrate,  and 
may  be  se;>arated  from  it  by  a buffer  layer. 

In  normal  FET  operation  the  source  and  drain  contacts  are 
connected  by  an  external  circuit  and  primary  current  flows 
through  the  conducting  channel.  The  width  the  depleted 
region,  and  consequently  that  of  the  conducting  chaivrvel,  is 
determined  by  the  potential  within  the  channel.  The  j)otentlal 
within  the  channel  Is  dependent  on  both  the  gate  and  drain 
potential  and  the  dc  resistance  is,  therefore,  bias  dei'endent. 
epical  source-drain  current-voltage  characteristics  are  there- 
fore nonlinear  and  dependent  on  the  gate  jmtential  (1). 

With  the  exception  of  very  recent  work  most  analyses  of 
FFT  devices,  e.g.,  amplifiers,  are  small  signal  studies.  The 
starting  point  for  these  studies  is  the  kmwn  dc  distribution 
of  space  charge,  electric  field  and  voltage  within  the  device. 
Theoretical  studies  usually  involve  twr -dimensional  approxima- 
tions to  the  three-dimensional  struetiire,  and  can  be  grf'uped 
into  either  analytical  or  numerical  approaches.  The  analyti- 
cal studies,  which  include  the  original  proposal  of  Ghnckley 
(1),  generally  consist  of  reducing  the  cr>m}>llcated  two-dimen- 
sional analysis  to  a series  of  separate  but  connected  one- 
dlmtnsional  problems.  The  analysis  I'eniiits  calculations  of 


the  current-voltage  characteristics,  the  transconductaace, 
etc.,  from  which  parameters  for  device  design  can  be  extracted. 

'nie  postulated  space  charge  distribution  within  an  FFT 
determines  hew  the  device  is  designed  for  a specific  applica- 
tion. The  adequacy  of  the  analytical  studies  for  determining 
the  space  charge  distribution  has  been  scrutinized  carefully 
as  inconsistencies  between  theory  and  experiment  have  emerged. 
In  one  general  study,  using  numerical  techniques,  Kennedy 
et  al,  (2)  predicted  that,  at  sufficiently  high  bias  values, 
regions  of  carrier  accumulation  and  depletion  would  form  with- 
in the  conducting  channel  of  the  device;  a phenomenon  not  pre- 
dicted by  earlier  analytical  studies.  The  presence  of  these 
nonuniform  space  charge  layers  was  shown  to  affect  the  gate 
bias  dependence  of  the  primary  current-voltage  characteristlca. 

Within  the  past  few  years  considerable  activity  has  cen- 
tered around  the  semiconductors  gallium  arsenide  Indium 
phosphide  as  FET  materials.  While  InP  is  of  recent  Interest, 
GaAs  FET's  provide  the  prospect  of  linear  amplification  of 
microwave  signals  of  up  to  at  least  10  GHz,  output  power  of 
several  watts  and  a bandwidth  of  an  octave  or  more.  From  the 
point  of  view  of  device  design  it  is  recognized  that  the  space 
charge  distribution  within  this  material  is  nonuniform  sind  in 
some  cases  electrically  unstable  (3).  The  reason  for  the 
latter  is  that  for  long  conducting  chamnels  and  sufficiently 
high  fields  the  electrons  may  acquire  enough  energy  to  become 
*hot*  and  create  a region  of  negative  differential  mobility 
where  the  average  mobility  decreases  with  increasing  electric 
field,  see  Fig,  2,  If  the  latter  occurs  the  device  becomes 
electrically  active  and  introduces  a mechanism  that  con^^etes 
with  the  small  signal  operation  of  the  device;  namely,  the 
spontaneous  nucleation  and  subsequent  propagation  of  local  re- 
gions of  high  electric  fields  (Gunn  domains). 

The  prospect  of  dealing  with  electrically  active  semicon- 
ductors introduces  additional  constraints  in  FET  design.  In 
the  first  place  low  resistance  contacts  are  known  to  allow  the 
formation  of  Gunn  domains.  In  certain  cases  instabilities 
associated  with  these  domains  can  be  eliminated.  In  other 
cases,  as  in  logic  applications,  Gunn  domain  nucleation  and 
propagation  are  features  incorporated  into  the  design  of  the 
device.  The  time  dependent  behavior  of  electrlcstlly  active 
semiconductors  cannot  always  be  described  by  small  signal  time 
dependent  analyses.  The  nature  of  the  problem  dictates  that 
the  transient  phenomena  be  treated  in  as  realistic  a way  as  is 


Vbd 


FUi.  1.  An  FCT  conflgur.tion  with  source  uid  drain  contacts 
on  the  ends  of  the  device  and  the  gate  contact  on  top.  Tlie 
device  is  modeled  in  the  xy  plane. 


FIELD  (VOLTS/CM) 

Fig.  2.  The  velocity  versus  electric- field  v(E)  curve  for 
n-G&As.  The  theoretical  curve  is  that  of  Butcher-Fawcett 
(6)  and  the  experiisental  points  were  obtained  by  Ruch- 
Kino  (7). 


for  carriers  of  charge  -e  where  e = l,6xlO“^9  coul.  < Is  the 
permittivity  of  the  material  which,  for  GaAs  is  9,8xlO“3  coul^/ 
Joule-cm.  V{x,y,t)  is  the  potential,  ?(x,y,t)  the  electric 
field,  N(x,y,t)  the  mobile  carrier  density,  and  I^(x,y,t)  the 
position  dependent  background  doping  density.  For  the  sasple 
calculations  NIq  is  a constant  equal  to  10^^/cm^. 

The  average  electron  velocity  is  related  to  the  electric 
field  through  the  relation 


numerically  possible.  Accurate  modeling  of  ',aAs  FPr*s  demands  w(x,y,t)  * - M ( |E  j )E(x,y,t)  (3) 

large  signal  time  dependent  analysis  in  at  least  two  dimen- 
sions. 


Large  signal  transient  analysis  of  FET's  has  been  done 
for  time  independent  potential  on  the  contacts  (4,5).  An 
initial  device  state  is  chosen  which  is  incoir^atible  with  the 
final  steady  state  and  the  resulting  transient  relaxation  is 
then  studied.  A more  difficult  but  more  realistic  problem  is 
solved  in  this  paper  in  idilch  the  semiconductor  equations  and 
the  external  circuit  equations  are  solved  «uid  coupled  together 
through  time  dependent  boundary  conditions,  specifically  the 
contact  potentials. 

The  principal  approximations  in  the  mathematical  formula- 
tion are  (a)  the  three-dimensional  carrier  flow  is  reduced  to 
two  dimensions  by  assuming  the  variables  to  be  independent  of 
the  z coordinate,  and  (b)  all  minority  conduction  is  ignored. 
Within  this  framework  a numerical  simulation  of  field  depen- 
dent nobility  devldes  is  performed.  The  mathematical  model  of 
the  semiconductor  consists  of  a set  of  nonlinear  partial  dif- 
ferential equations  which  are  solved  on  a rectangular  mesh 
using  explicit  finite  difference  techniques.  Particular 
attention  is  paid  to  the  selection  of  boundary  conditions. 

The  boundary  conditions  approximate  the  physical  and  electri- 
cal characteristics  of  the  outer  periphery  of  the  semiconduc- 
tor structure.  The  exposed  surfaces  of  the  semiconductor  are 
assumed  to  be  ideal  electrical  insulators;  and  no  current  is 
permitted  normal  to  these  bounding  surfaces.  The  low  resis- 
tance contacts  and  the  gate  contact  are  approximated  by 
equlpotential  surfaces  with  a prespecified  charge  density. 

For  the  calculations  illustrated,  the  low  resistance  source 
€Ukl  drain  contacts  are  charge  neutral  and  located  sufficiently 
far  from  the  active  region  of  the  device  so  that  for  the  cur- 
rent levels  involved  they  have  no  influence  on  the  electrical 
properties  of  the  device,  A specific  impurity  atom  distribu- 
tion is  assumed,  which  in  general  can  be  spatially  dependent. 
In  the  sample  calculations  it  is  assumed  to  be  homogeneous. 

The  semiconductor  equations  are  coupled  to  the  circuit 
equations  by  the  boundary  conditions  at  the  device  terminals. 
The  solution  of  the  resulting  time  dependent  boundary  value 
problem  is  a feature  of  the  simulation  and  rejvesents  an 
advance  in  the  state-of-the-art  of  these  types  of  calcula- 
tions. 

The  simulation  is  implemented  using  interactive  conqxiter 
graphics.  As  the  solution  evolves,  displays  of  the  potential, 
charge  density  and  current  flow  in  the  semiconductor  are 
available  via  a computer  graphics  terminal. 


where  w(x,y,t)  is  the  average  electron  velocity  and  m(  E ) is 
the  field  dependent  mobility. 

The  equation  for  carrier  transport  including  diffusion  is 
T(x.y,t)  = eN(x,y,t)n(lEI)E(x,y,t)  * eDgr^l  N(x,y,t)  (U) 


where  D is  the  diffusion  coefficient,  which  is  generally  elec- 
tric field  dependent.  In  the  saa^le  calculations,  D is  taken 
to  be  constant  and  equal  to  200  cm^/sec.  above  equations 
are  subject  to  the  constraints  of  the  equation  of  continuity 

div  J(x,y,t)  = e [N(x,y,t)  - IlQ(x,y,t)|  (5) 
or 


dlv  T(x,y,t)  » 0 (6) 

where 


-T(x,y,t)  = J(x,y,t)  ♦ (7) 


The  boundary  conditions  to  the  above  equations  are  de- 
rived from  the  surface  and  contact  properties:  no  current 
flows  through  the  free  surfaces.  With  the  assuinjtlons  of 
zero  permittivity  for  the  surrounding  space,  this  leads  to 

ngrad  V > ngrad  N > 0 (8) 

along  the  free  surfaces,  n is  a unit  vector  normal  to  the 
free  surface  of  the  semiconductor.  At  the  source,  drain  and 
gate  contact  N(x,y,t)  is  preaaaigned.  For  the  saiiq>le  calcula- 
tions 


farnula 


N = , 

source  drain 


Nq  = 10^^/ca^ 


(9) 


and 


V.  +V  +V  +V  -Uv 

k.l+l.m  k,l-l,in  k,t,in+l  k,l,itt.l  k.l.s 

= Vi.»’ 


(16) 


gate 


= iO"^  N 


(10) 
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The  source  potential  is  set  to  zero  and  the  gate  potential  Vq 
and  drain  potential  Vjj  are  determined  by  simultaneous  solution 
of  the  semiconductor  equations  and  the  circuit  equations. 

For  the  three  contact  device  of  Fig.  1 the  circuit  equa- 
tions are 


(U) 


The  boundary  potentials  V^.I^l.rn.  ^k.i.o,  Vj,  , are 

determined  from  Dirichlet  conditions  on  the  contacts  and 
Neumann  conditions  elsewhere.  The  potential  on  the  contact  is 
found  by  solution  of  the  external  circuit  equations.  Away 
from  the  cont^t  the  normal  electric  field  is  assuoed  to  be 
zero  on  the  boundary. 

The  con^nents  of  electric  field  are  defined  between 
nodal  points  as  follows: 


(12) 

* yt) 

(13) 

l3(t)  = I^(t)  * 

yt) 

(lU) 

*k, I ,m 


(V  - V ) 

^ ky  i,m  k,,|wl,m^ 


l,m 


I,m  ' l,m-l^ 


(IT) 


where  Vs  is  the  source  potential,  Vq  the  gate  potential,  Vp 
the  drain  potential,  I5  the  source  current,  Iq  the  gate  cur- 
rent, If)  the  drain  current,  Zq  the  gate  impedance,  Zj,  the 
drain  impedance.  The  external  circuit  impedances  are  modeled 
as  liiBq)ed  linear  elements.  The  current  passing  across  a coa> 
tact  is  the  integral  over  the  contact  area  of  the  component 
of  current  density  normal  to  the  contact.  For  a device  of 
dimension,  1^,  in  the  z direction,  euid  a contact  extending 
over  a region,  f , 

I(t)  - lij^  n T(x,y,t)  d('  (15) 

The  loop  currents  and  external  circuit  potentials  are  normal- 
ized with  respect  to  the  peak  average  electron  velocity  for 
"/aAs  and  the  corresFonding  electric  field,  respectively.  The 
loop  cxirrents  and  circuit  impedances  are  further  normalized 
with  respect  to  the  device  depth, 


As  defined  above  I is  located  midway  between  the  nodes 
(i-l,m)  and  wfille  Eyj^  I is  located  midway  between 

nodes  ( i,m-l)  and  (i,m).  In  order  to  coo^te  the  mobility, 

E ),  an  average  electric  field  E(av)  calculated  whose 
X and  y components  are 


^(»v)k,  |,B  “ °-5 

(Ifi) 

The  vector  E(av)  computed  is  located  in  the  center  of  a 

square  which  has  at  its  corners  the  nodes  (1,111),  (l-l,m), 

( i,nk-l),  (i-l,ro-l).  Numerical  experiments  have  indicated 
that  the  method  of  computation  of  E(av)  s^Tects  the  stability 
of  the  solution.  Other  methods  which  l».ve  been  reported  (U,5) 
lead  to  more  conplex  or  less  stable  algorithms  than  Eq.  (18). 
The  average  electron  velocity  is  confuted  from  ^(av)  ^ 
follows : 


The  sign  of  the  current  density,  i,  as  well  as  the  signs 
of  the  loop  currents  in  the  external  circuit  are  consistent 
with  conventional  current  flow.  In  particular,  for  the  cir- 
cuit of  Fig.  1,  positive  source  current,  I5,  implies  electron 
flow  into  the  device  at  the  source  contact  and  positive  drain 
current.  Ip,  implies  electron  fl<-v  out  of  the  device  at  the 
drain  contact.  Positive  gate  current,  Iq,  Imxilles  electron 
flcar  out  of  the  device  at  the  gate  contact. 

DIFFtyewCE  EQUATIONS 


The  device  potential  and  mf>bile  carrier  density  are  conw 
puted  on  a uniform  rectangular  mesh.  The  mesh  contains  LxM 
internal  nodes  with  bf^undary  conclitlons  lmi-^)sed  via  addi- 
tional nodes  at  the  edges  of  the  device.  Therefore,  the 
total  number  of  nodes,  Including  the  boiindTiry,  is  (h*2)x(M+2). 
Derived  quantities  such  as  electric  flej  i,  chfirgc  ontl  t<;tal 
current  densities  are  defined  at  j*olnt:  other  than  the  nodes 
although  they  are  written  In  terms  •>1'  the  ryxifj  indices.  Let 
I be  the  index  in  the  y dlrectlim,  m the  Index  in  the  x direc- 
tion, and  X the  uniform  mesh  increment,  lot  K !>e  u time 
Index  and  T the  time  Increment. 

frisson’s  equation  is  dlrcrctlzed  usirv'  the  five  jriint 


’'k.l.m"  *''^l®(«v)l^^(»v)k,  l,B 


(19) 


BO  that  7 la  defined  at  the  seine  location  as  E(av).  k spe- 
cific veloeity.electric  field  curve  is  included  in  the  simu- 
lation via  a piecewise  linear  map  of  nobility  versus  electric 
field  magnitude.  The  carrier  tranaport  equation  is 


^k,l,m  ^K,l,m+1^ 


CD 


X I ,m  I -l,m^ 


1 =1  I s I/+1 
IsmsM 
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^ 1,111  * l,ii>-l^^*yk,  I.m  * *yk,i*l,B^ 
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(20) 


eD 


X ^**k,  l,m  " l,iD-l^ 


Is  isL 
1 s n s M+1 


The  charge  current  density,  J,  Is  defined  such  that  Jxk  | q 
is  located  midway  between  nodes  ( l,m)  and  (l-l,m)  and 


>Iyk,  1,111  located  nidway  between  nodes  ( l.m)  and  ( l,»-l). 

An  explicit  integration  fonula  is  used  to  propagate  the 
snblle  carrier  density  in  time- 


nie  solution  requires  initial  conditions  >o  I in'  ^ > ^Do 
and  boundary  conditions  on  V and  N.  At  each  tine  step  tSe 
solution  procedure  may  be  suanarlzed  as  follows,  given  | g,, 
^k-l*  ^Dk-1  estimates  of  Iqj^,  ’ ’ 


*Sc+l,  I ,m 


'V,l,m  * ex  *^*k,Ul,m  ' i,m  * '^yk,  I .m-l  ‘ '^yk,l,m^ 
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(21)  Relaxation  iteration  is  used  to  update  the  estijnates  of  the 

total  currents  I(^  and  When  the  iterations  have  converged 

the  mobile  charge  density  is  propagated  to  the  next  time  step 
using  Eq.  (24). 


where  Nq  has  been  assumed  to  be  tijne  independent. 

The  total  current  density,  T,  is  calculated  so  that  its 
components  are  defined  at  the  same  locations  as  the  charge 
current  density: 

T ^^k+l,l,ni  ^k,l,ni^ 

1 - I ■ L+1 
1 ' m * M 

(22) 

^k,  I ,m  " *'^yk,  I ,m  ~ T ^^k+1, 1 ,in  " ^k,  I ,in^ 

15  15  L 
1 5 m 5 M+1 

The  total  current  on  each  contact  is  founvi  by  integrating  the 
normal  component  of  total  current  density  over  al ^ nodes 
associated  with  the  contact. 


The  solution  of  the  difference  equations  is  straightfor- 
ward with  the  exception  of  the  discrete  farm  of  Poisson's 
equation,  Eq.  (l6).  Methods  for  the  solution  of  Poisson's 
equation  are  summarized  and  compared  by  Dorr  (8).  From  the 
available  techniques  direct  solution  by  Gaussian  elimination 
has  been  chosen.  Gaussian  elimination  is  not  generally  used 
when  two  space  dimensions  are  involved  since  precos^utation 
and  storage  of  an  LxbcM  array  is  required.  However,  when  the 
solution  mesh  is  such  that  L is  much  less  than  M,  as  is  the 
case  with  the  devices  being  studied,  and  when  the  mesh  size  is 
not  too  large,  Gaussian  elimination,  by  virtue  of  its  con5)uta- 
tionaJ.  efficiency  for  time  dej'endent  problems,  is  the  method 
of  choice.  The  procedure  for  solution  of  the  discrete 
Poisson's  equation  in  two  dimensions  by  'Gaussian  elimination 
is  reviewed  in  (8). 

Stability 

Stability  analysis  of  the  device  difference  equations  is 
j'ossible  only  when  linearizing  assumptions  are  made.  * The  sta- 
bility of  the  explicit  integration  of  Eq.  (21)  can  be  shown  by 
von  Neumann  stability  analysis  to  be  limited  by 


The  external  circuit  equations  are  solved  using  an 
iji^liclt  scheme.  For  resistive  circuit  elements 


% " ''BDk  - Bd  • Ir>k 
Vo,  = - R,  . 


(231 


where  Rjj  is  the  irain  circuit  resistance  and  R,.  is  the  gate 
circuit  resistance. 

Method  of  oclution 

The  presence  of  the  external  circuit  intr^vtuces  some 
difficulties  into  the  solution  of  the  equations  describing 
the  device  charge  distribution.  The  difference  equations 
have  been  derived  assuming  time  indejendent  boundary  poten- 
tial. Referring  to  Eq,  (22),  unless  the  boundary  potentiaLL 
at  time  k*l  is  the  same  as  the  boundary  j'otential  at  time  k, 
the  current  density  will  not  be  correctly  ct)mi>uted.  In 
particular,  the  current  density  in  the  region  of  the  contacts 
will  not  satisfy  Eq.  (6);  that  is,  it  will  not  be  divergence 
free.  The  difficulty  is  overcome  by  the  use  of  the  foUt-wing 
artifice:  the  mobile  carrier  density  is  pro}.»agated  forward 
in  time  by  a small  increment,  T’  T,  •’ver  which  it  is  suffi- 
ciently accurate  to  assume  that  the  boundary  potential  does 
not  change.  Let  be  the  mobile  carrier  density  which  has 

been  propagated  by  T’  • Then  ]iro{»agatlon  to  the  next  full 
time  step  is  performed  as  follows: 


T 

**k*l,  I ,m  ^k,  I ,m  * ^”k*ft,l,m  \,l,m^ 


(24) 
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4d 


(25) 


when  constant  boundary  potential  and  constant  average  electron 
velocity  ( vT  - -v  = constant)  are  ass\imed,  ^rplcal  values  for 
JaAs  are  X = 10"^  cm,  D = 200  cm^/sec,  and  v = 10^  cm/sec  so 
that  the  stability  limit  is  T x2/4D  = 1,25x10“^^  sec.  Note 
that  for  l(>w  values  of  diffusion,  stability  is  limited  by 
T D/v^.  Implicit  integration  of  carrier  density  has  been 
used  in  simulations  of  silicon  semiconductors  which  exhibit 
IcTw  diffusion  coefficients  (4).  Numerical  experiments  have 
indicated,  however,  that  for  dc'/ices  with  field  dependent 
mobility,  #i{  E ),  the  stability  criteria  derived  from  the 
linearized  system  are  of  limited  usefulness  in  predicting  per- 
f<3rmance  of  a particular  method.  The  sample  calculations  to 
be  presented  subsequently  were  performed  at  a time  increment 
of  3»9^>x10"^3  sec,  numerical  stability  was  maintained  for 
the  sane  time  increment  when  the  diffusion  coefficient  was 
reduced  to  zero.  Additional  numerical  experiments  have  shewn 
that,  for  laAc  devices  operating  in  the  negative  differential 
mobility  regime,  implicit  integration  affords  no  iitqjrovenent 
in  computational  efficiency  over  explicit  integration. 

Analysis  of  the  stability  of  the  external  circuit  equa- 
tions is  pi^ssible  when  the  eleotric  field  in  the  device  is 
assumed  to  be  homogeneous  and  of  sufficiently  low  magnitude 
that  average  electron  velocity  is  a linear  function  of  elec- 
tric field.  When  these  assumptions  are  made  the  device  is 
seen  to  be  capacitive.  Thus,  when  the  external  circuit  Iji^ed- 
ances  are  resistors  a one  time  step  delay  applied  to  the  loop 
currents  in  the  external  circuit  equations,  which  leads  to  an 
explicit  exiiression  for  the  boundary  potentials,  also  leads  to 
numerical  Instability.  Therefore,  the  iji^jllcit  fonmilas, 

Eq.  (23),  are  emi'loyed.  A relaxation  iteration  with  a relaxa- 
tion constant  less  than  t>ne  (under  relaxation)  is  used  for 
both  the  gate  and  strain  circuits  (9). 


SAMPLE  CAI£UIATIOKS 


The  slmlatlon  Is  demonstrated  by  tvo  sample  calculations 
presented  in  Figs.  3-7.  The  solution  Is  found  on  an  8x40  rmsh 
exclusive  of  boundary  nodes.  The  configuration  which  Is  simu- 
lated Is  that  of  Fig.  1 In  which  the  source  and  drain  are  on 
the  sides  of  the  device  and  the  gate  Is  on  the  top.  The  param- 
eters are  as  follows: 


device  length  (x) 
device  width  (y) 
X 
T 


gate  width 


2d 

D 


10  microns 
2.19  microns 
0.244  microns 
3.96x10"^3  sec 
3.416  microns 
0 

Ro  . 

200  cm/sec*^ 
10^5  cm" 3 


by  ap})lying  a small  drain  bias  potential  to  the  device  and  run- 
ning the  sjjnulatlon  until  the  mobile  carrier  density,  irhich  Is 
initially  set  to  reaches  steady  state. 

Figure  3 shows  the  time  and  space  evolution  of  the  mobile 
carrier  density  and  the  vector  distribution  of  current  density 
within  the  FCT  which  occur  when  the  drain  bias,  V5]),  Is  raa9>ed 
from  0.1  to  1.0,  in  normalized  units ^ in  4 picoseconds,  ^e 
left-hand  column  shows  N(x,y,t)  at  successive  instants  of  time. 
For  clarity  in  the  display  N(x,y,t),  expressed  in  units  of  Hq, 
Increases  in  the  downward  direction.  The  current  density  dis- 
plays are  vector  displays  shewing  both  magnitude  and  direction. 
The  line  segments  reduce  to  points  when  the  magnitude  of  the 
total  current  density  is  zero  and  reach  their  maxiiauiti  length 
when  the  magnitude  of  the  current  density  is  greater  than  or 
equal  to  the  j^eaX  current  for  GaAs  (set  Fig.  2).  The  time 
evolution  of  the  loop  currents  and  the  drain  potential  are  dis- 
played in  Fig.  4 in  normalized  units. 


where  Rq  is  the  intrinsic  low  field  resistance  which  depends  on 
the  dimensions  of  the  device  in  the  z direction.  Given  the 
loop  current  normalization  with  respect  to  device  depth,  It  is 
sufficient  to  state  that  Zq  = Rq.  The  background  charge  den- 
sity on  the  gate  conteict  Is  set  to  10“ Nq.  A three  piece 
linear  approximation  of  the  velocity-electric  field  curve  of 
an  NPM  device  has  been  used  in  which  the  electron  velocity  at 
high  field  is  one  half  the  peaX  velocity.  The  source  and  gate 
potentials  are  zero.  Realistic  Initial  conditions  are  achieved 


The  details  of  the  disple^  show  an  evolving  depletion 
layer  that  moves  into  the  semiconductor  under  the  gate  elec- 
trode. The  depletion  layer,  where  N{x,y,t)<NQ,  iji5)lles  a 
tranrverse  comjKjnent  of  electric  field  pointing  in  the  direc- 
tion of  the  gate  electrode  and  causing  electron  flow  into  the 
r.emicnnductor  through  the  gate.  Note  the  current  density  dls- 
;1^  of  Fig.  3 at  ^.3  ps  and  also  Fig.  4 where  the  magnitude 
of  the  gate  current  initially  exceeds  the  source  current.  In 


Fig.  3.  The  time  and  space  evolution  of  the  normalized  mobile 
carrier  density,  N/I^,  and  ciirrent  density,  i,  which  t^ccurs 
when  the  drain  bias  potential  is  ram|»ed  from  0,1  to  1.0  in 
4 pe. 


Fig,  4.  The  time  evolution  of  the  loop  currents  and  drain 
potential  when  the  drain  bias  potential  is  raugied  from 
0.1  to  1.0  in  1*  ps. 


steady  state  the  gate  current  approaches  zero.  It  is  signifi- 
cant to  note  in  Fig.  3 that  the  current  density  in  steady  state 
is  greatest  under  the  gate  contact  where  the  cross-sectional 
area  of  the  conducting  channel  is  least.  For  the  particular 
set  of  bios  ]x)tentials  all  field  values  within  the  conducting 
channel  are  less  than  the  electric  field  at  peaJc  current  for 
laAs  and  therefore  the  increase  in  current  density  is  accoirw 
plished  without  any  local  regions  of  charge  accumulation. 

In  the  second  sequence,  shi-wn  in  Figs.  5-7,  the  drain 
bias  jK'tential  is  rain)«d  from  0.1  to  2.0  while  the  source  and 
gate  potentials  are  held  at  zero.  In  the  absence  of  a gate 
electrode,  such  a drain  bias  potential  would  not  be  suffi- 
ciently large  to  cause  the  device  to  operate  in  the  region  of 
negative  differential  mobility.  fKvever,  the  constriction 
caused  by  the  gate  electrode  increases  the  magnitude  of  the 
electric  field  under  the  gate  to  a level  such  that  the  device 
oi^erates  locally  in  the  NDM  regime.  The  scale  of  Fig.  5 is 
such  that  the  length  of  tJje  line  segments  rej^resenting  ciirrent 
density  merge  into  a s»’lid  line  just  as  the  NDM  regime  is 
entered.  The  area  in  which  NDM  operation  is  initiated  may  be 
clearly  seen  just  unler  the  gate.  An  accumulati^'n  of  mobile 
carriers  occurs  and,  in  the  sequence  sli'wn  in  Fig.  6,  nv'ves 
from  the  gate  to  the  strain  where  the  accumulation  <iecays.  The 
process  then  repeats  itself  an^i  will  <U‘  so  indefinitely.  The 
evidence  of  oscillations  within  the  ievice  is  the  presence  of 
oscillations  in  the  l<K'p  current  in  the  external  circuit. 

Fig.  7.  lii'te  in  Fig.  7 that  the  drain  potential  is  time  de- 
pendent and  therefore  an  ass*ijnption  of  ct'nstant  piotential 
would  not  be  realistic  f‘’r  this  situatiem. 


N/No 


A c*'mputer  simulation  has  been  presented  which  all«ws 
large  signal  tine  dependent  analysis  cf  b’jIX  semiconductors 
oierating  in  an  external  circuit.  Hie  simulation  was  develori- 
ed  to  stuij,'  negative  differential  mobility  devices;  h'vever, 
any’  semiconducting  material  for  which  a velocity-electric 
field  curve  is  known  may  be  modeled.  Alth< ugh  the  sample  cal- 
culations presented  were  p'erfermed  f » r a sp>ecific  contact  con- 
figuration, the  simulation  p-ermits  placement  of  the  contacts 
anywhere  on  the  boundaries.  The  primary  limitation  of  the 
numerical  technique  is  that  the  device  must  be  modeled  rn  a 
rectangle. 

By  simulating  the  interaction  <if  the  semiconducting  de- 
vice and  the  external  circuit  a m re  representative  analysis 
of  ievice  behavior  is  possible;  particularly  in  those  situa- 
tions where  the  assupq'tion  of  Constant  electrode  potential  is 
not  valid. 

This  study  was  sup.p’O'rted  by  the  . Ilavy,  •'ffice  of 
Naval  Research, 


N/Nq  current  density,  i 


Fig,  5.  The  normalized  mr-bilc  carrier  density  and  c^irrent 
iensity  at  the  initiation  • f an  accumu.lati'»n  of  mobile 
carriers 


Fig.  (■>  A time  sequence  showing  the  formation  and  prop'au’a- 
tion  of  dunn  domains  in  a PaAs  FKT 
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Fig.  7.  The  time  evolution  of  loop  currents  and  drain  poten- 
tial which  occurs  during  Gunn  type  oscillations.  The  tran- 
sient  euid  the  first  cycle  of  the  emcillation  is  shown. 
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